Vortex-line solitons in a periodically modulated Bose gas 
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We study the nonlinear excitations of a vortex line in a Bose-Einstein condensate trapped in a 
one-dimensional optical lattice. We find that the classical Euler dynamics of the vortex results in 
a description of the vortex line in terms of a (discrete) one-dimensional Gross-Pitaevskii equation, 
which allows for both bright and gray soliton solutions. We discuss these solutions in detail and 
predict that it is possible to create vortex-line solitons with current experimental capabilities. 



Introduction. — Solitons are special solutions of non- 
linear differential equations with particle-like properties. 
In a Bose-Einstein condensate with a repulsive interac- 
tion a gray soliton corresponds to a density minimum for 
which the loss in interaction energy is balanced by the 
increase in the kinetic energy. Gray solitons in a Bose- 
Einstein condensate have been observed experimentally 
by imprinting a ir phase step onto the condensate wave 
function 0, Q ■ While a Bose-Einstein condensate with 
a repulsive interaction is stable against collapse, a con- 
densate with attractive interactions may collapse. This 
collapse in first instance results in ever increasing den- 
sity gradients. However, at later stages of the collapse 
the attractive interaction energy may be balanced by the 
kinetic energy, which tends to spread out the condensate 
wave function. Therefore, in a one-dimensional conden- 
sate a collapse leads to a creation of a train of bright soli- 
tons. Bright solitons in a Bose-Einstein condensate have 
indeed been observed experimentally 0,0]- Also theoret- 
ically the collapse of the Bose-Einstein condensate into a 
train of bright solitons has been studied 0, . 

In this Letter we point out the possibility of vortex- 
line solitons in a Bose-Einstein condensate in a one- 
dimensional optical lattice. As is well known, the physics 
of a dilute Bose gas in an optical lattice is described 
by the Bose- Hubbard model 0, @- which in particu- 
lar predicts a quantum phase transition from a super- 
fluid state into a Mott-insulator state when the lattice 
becomes sufficiently deep. This transition was recently 
observed by Greiner et al. || . In the situation of inter- 
est to us here the one-dimensional optical lattice splits 
the Bose-Einstein condensate into a stack of weakly- 
interacting pancake condensates, and the vortex line in- 
tersects each one of them. We show that due to the 
inhomogeneous density distribution of the Bose-Einstein 
condensate, the dynamics of the vortex line is governed 
by a (discrete) one-dimensional nonlinear Schrodinger or 
Gross-Pitaevskii equation with the familiar soliton solu- 
tions. 

These solitons are possible when the pressure due to 
the kinetic energy balances the pressure due to the in- 
teraction energy. In principle, these competing pro- 
cesses also exist for a vortex line in an ordinary magneti- 



cally trapped Bose-Einstein condensate without an opti- 
cal lattice. In particular, the inhomogeneous condensate- 
density distribution in the radial direction gives rise to 
an interaction-like energy term and the energetic cost for 
bending the vortex line gives rise to the analog of the ki- 
netic energy. However, in an ordinary Bose-Einstein con- 
densate the vortex stiffness is orders of magnitude larger 
than in an optical lattice, where the vortex stiffness can 
be reduced by simply increasing the lattice depth. In 
practice this implies that the size of the soliton in an or- 
dinary condensate is very large and, in particular, larger 
than the typical condensate size. Setting aside other com- 
plications this fact alone suggests that with current ex- 
perimental capabilities vortex-line solitons are a possibil- 
ity that only exists for a Bose-Einstein condensate in an 
optical lattice. 

So-called optical vortex solitons llfj and solitonic vor- 
tices in a Bose-Einstein condensate fjjj have been dis- 
cussed before. Neither one of them is, however, related 
to the vortex-line solitons we discuss here. The former 
are related to the non-spreading propagation of a light 
beam with a vortex phase pattern through a nonlinear 
medium, while the latter refers to the soliton-like prop- 
erties of a vortex moving through a very narrow chan- 
nel. Moreover, in classical fluid dynamics Hashimoto 0] 
mapped the vortex line dynamics, in the so-called local 
induction approximation, into a nonlinear Schrodinger 
equation. As a result, solitons predicted by this formu- 
lation are only mathematically related to our vortex-line 
solitons and the physics is fundamentally different. 

Theory. — We assume a Bose-Einstein condensate ex- 
periencing a one-dimensional optical lattice potential in 
the longitudinal direction. This lattice potential splits 
the Bose-Einstein condensate in N s sites and each site 
has N atoms. In the radial direction the atoms are mag- 
netically trapped by a harmonic trap with a trapping 
frequency ui r . The magnetic trapping in the longitudi- 
nal direction is assumed to be so weak that it can be 
safely neglected. While the lattice is taken to be deep 
enough to allow us to use a tight-binding approximation 
and to include only the weak nearest-neighbor Josephson 
coupling, it is also taken to be shallow enough to sup- 
port a superfluid state as opposed to the Mott-insulator 
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state |3. 

In our earlier work [Til we showed that up to sec- 
ond order in the vortex displacements from the center 
of the Bose-Einstein condensate the Hamiltonian for the 
vortex line is that of a set of coupled harmonic oscil- 
lators. The eigenmodes of this system are the Kelvin 
modes. These modes have been recently observed in a 
cigar shaped Bose-Einstein condensate in the absence of 
an optical lattice [lfi, Il7j . However, if the vortex en- 
ergy functional is expanded up to fourth order in the 
displacements, the physics of the vortex line turns out to 
be described by a one-dimensional Bose-Hubbard model. 
The corresponding Hamiltonian is given by 



H = £^ [fiwo + Jv] fi+t 
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where (n, m) denotes the nearest-neighbor layers and the 
Kelvon operators v n were defined in terms of the vortex 
displacements (x n ,y n ) {R/2VN)(vl + v n ) and 

y n = {iR/2yN){v\ — v n ), where R is the radial size of 
the pancake condensate. 

In view of the complicated dynamics of the three- 
dimensional vortex line in a bulk superfluid, this rela- 
tively simple result is quite remarkable. Furthermore, as 
mentioned previously the Bose-Hubbard model shows a 
quantum phase transition from a superfluid into a Mott- 
insulator state when Jy /Vq becomes sufficiently small. In 
our case this transition describes the break-up of the vor- 
tex line into individual pancake vortices due to quantum 
fluctuations. However, in the following we work solely 
in te superfluid regime, where a classical mean-field de- 
scription is appropriate. 

Using a Gaussian variational ansatz for the condensate 
wave function in each site j3| that is proportional to 
exp [— (x 2 + y 2 ) /2i? 2 ] , we can analytically calculate the 
strength of the nearest-neighbor coupling to be 
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and the interaction strength 

2huj r (l r /R) 2 T [0, Z 4 /i? 4 ] - 3hcj r (l r /R) 2 - AhVt 



Vn = 



47V 



'8tt 2 Vl 1 'm\ 2 is the oscillator frequency 
r [a, z] is the in- 



where lul 

of the optical lattice, l r — ^Jh/mui r , 
complete gamma function, Vl cos 2 (2ttz/\) is the lattice 
potential, A is the wavelength of the laser light, and f2 
is the rotation frequency of the magnetic trap. Fur- 
thermore, a straight and slightly displaced vortex pro- 
cesses around the condensate center with the frequency 
loq = (io r l 2 r /2R 2 ) (1 -T [0,^ 4 /i? 4 ]) + Q. The radial size 



R of the Bose-Einstein condensate is determined by min- 
imizing the condensate e nerg y with the vortex line in the 
center of the condensate |l5j . 

In Ref. [? ] we discussed the equilibrium properties of 
a straight vortex when the interaction term is attractive 
and pointed out that in this case the quantum mechan- 
ical state of the vortex can become strongly squeezed. 
However, here make the usual classical approximation 
and study the nonlinear excitations of the vortex line. 
The solitonic solutions that we find can in principle also 
be squeezed under the appropriate conditions, but this 
interesting topic is deferred to future work. 

Under the assumption of a coherent state for the Bose- 
Einstein condensate, the Euler dynamics of the vortex 
line in an optical lattice is given by the discrete Gross- 
Pitaevskii equation 



iv n (t) 



■h 



K-i(i) - 2u n (i) + v n+1 (t)} 



[uj + V \v n (t)\ 2 ] v n (t), 



(2) 



where from now on we use units such that % = 1. If the 
variation of the vortex positions is small over one lattice 
spacing this equation is simply the discretized version of 
the continuum Gross-Pitaevskii equation 



iv n (t) 



Jy d 2 v n {t) 
~2 dn 2 



+ [u Q + V \v n (t)\ 2 ]v n (t). (3) 



As indicated by the above equation we from now on also 
choose to work with a longitudinal unit of length that 
equals the lattice spacing A/ 2. 

Dark vortex-line soliton. — Dark or gray soliton solu- 
tions of Eq. © are possible when the interaction is re- 
pulsive, i.e., Vo > 0. This is the case in the non-rotating 
trap, for example. Explicitly the gray soliton solution 
moving with velocity v is given by 



v n {t) = \/no < i cos 9 + sin 9 tanh 



sin 



(n - vt) 



-ifit 



where uq = linin^oo |f n | 2 , cos6> = v/c, c = V JvVqUq 
is the speed of sound, fi = loq + Vq'^Oj an d the soliton 
size is £ = sj Jy/Von lattice spacings. While the gray 
soliton solution is simple to construct theoretically, its 
experimental preparation is challenging since it requires 
a good control over the position of the vortex in each 
site. One possible way to create a dark soliton is to 
prepare a vortex line that is tilted with respect to the 
optical potential. If the vortex line crosses the center of 
the condensate somewhere along the lattice, this initial 
state provides a prototype for the dark soliton. Other 
initial states can similarly resemble gray solitons. Un- 
fortunately, the ensuing time evolution, which we have 
studied numerically, typically also generates a fair num- 
ber of other vortex- line excitations, such as Kelvin modes 
and phonons. This makes the experimental observation 
complicated. 
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FIG. 1: A bright soliton solution of Eq. |QJ with N = 1000, 
k = 0, JV = 0.4 fiwr, and Vb = — 7-10 hui r . These values for 
Jv and Vb are realistic in an optical potential with a lattice 
spacing d = 398 nm, depth Vl = 15 E r , with 100 87 Rb atoms 
in each site, and a rotation frequency £1 = 0.2 lo t . here uj r = 
2iy ■ 100 Hz is the radial trapping frequency, E r = h 2 /2mA 2 
is the recoil e nergy, a nd the unit of length in the transverse 
direction is ^Jh/muj r . 



Bright vortex-line soliton. — Since experimental prepa- 
ration of the dark soliton is complicated, we focus now on 
the regime where the interaction is attractive, i.e, Vq < 0. 
This occurs, when the magnetic trap is rotated so fast 
that the vortex is energetically stable. With an attrac- 
tive interaction there exists a solution of Eq. (J3J) corre- 
sponding to the bright vortex-line soliton moving along 
the lattice, which is given by 



v n (t) 



2£ cosh[(n- J v kt/2) /£] ' 



(4) 



where k is the wavevector, \i = ui n — Jy/2£ 2 + Jyfc 2 /2, and 
the size of the soliton is £ = 2JV/|Vo|iV lattice spacings. 
This solution is normalized such that J dn\v n \ 2 = Nq. 
When k = 0, this solution is equivalent with a shape pre- 
serving precession of the curved vortex line. Since the 
chemical potential /j, of the stationary soliton is smaller 
than loq, the soliton precesses around the center of the 
condensate with a lower precession frequency than a 
straight vortex line. In Fig. ^ we show an example of 
the bright vortex soliton with typical parameter values. 

Soliton creation. — In order to gain quantitative under- 
standing of the vortex-line collapse into a train of such 
solitons, we need to study the excitations around the ho- 
mogeneous solution. At long wavelengths the Bogoliubov 
dispersion for a displaced but straight infinitely long vor- 
tex line corresponding to the density \v n \ 2 = uq is given 



At wavelengths that are comparable to the lattice spac- 
ing the lattice discreteness becomes apparent and the 
dispersion is modified by the replacement Jyk 2 /2 — > 
Jv [1 — cos(fc)]. However, for our purposes the long wave- 
length limit of the dispersion is sufficient. When Vb < 
the dispersion has imaginary energies and the system is 
dynamically unstable. The most unstable mode has the 
wavevector fco = y/2\Vo\n /Jv and the absolute value 
of its energy can be used to estimate the characteristic 
timescale r for the collapse by t = l/|efc | = 1/| VojnQ. 
In contrast to the nearest-neighbor coupling strength Jy , 
Vb is not very sensitive to the depth of the lattice poten- 
tial. Therefore, the timescale for the vortex-line collapse 
is hard to tune by changing the lattice depth. A better 
way to tune the collapse time is to vary the initial vor- 
tex displacement. We demonstrate this in Fig. [3 which 
shows the behavior of r with some typical parameter val- 
ues as a function of the initial vortex- line displacement. 
It is clear from this figure that by displacing the vortex 
line sufficiently, the collapse can occur so quickly that 
the dissipative processes are unlikely to affect the vortex 
dynamics. 

As the wavevector fco = 2-7r/Ao indicates the wave- 
length Ao of the most unstable vortex-line fluctuations, 
we can estimate that the collapse creates about N s ko/2TT 
solitons. This in turn can be used to estimate the size of 
the solitons after the collapse. We find, using the lattice 
spacing as the unit of length, that the soliton size after 
collapse is given by £ = ^/2Jy/|Vb|n /7r. However, it is 
clear that the above arguments are only valid if the fi- 
nal soliton size is much smaller than the system size N s . 
Otherwise, boundary effects are expected to influence the 
results considerably. Another way to understand this is 
to note that in the finite lattice the smallest nonzero wave 
vector is given by k m i n — 2tt/N s . Moreover, we can see 
from Eq. JSJ that the dynamical instability exist only in 
the regime where |fc| < 2y/\Va\no/ Jv = k c . Therefore, 
the vortex line in an optical lattice can only be dynam- 
ically unstable if k m i n <C k c . This implies the condition 




«1, 



(6) 



which turns out to be equivalent to the requirement that 
the final soliton size must be much smaller than the sys- 
tem size. 

The above arguments based on a Bogoliubov approach 
are valid only when the deviations from the initial ho- 
mogeneous system are small. Given enough time the dis- 
turbances in the dynamically unstable system grow large 
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FIG. 2: The timescale for the vortex line collapse as a function 
of the initial displacement ro relative to the condensate size R. 
For concreteness we further assumed that 1000 S7 Rb atoms 
were trapped in each site with Q — 0.2 u) r and ui r = 27r- 100 Hz. 
The lattice depth was Vl = 15 E r and the lattice spacing was 
d = 398 nm. 



and the Bogoliubov approach fails. In this limit we have 
to solve the discrete time-dependent Gross-Pitaevskii 
equation numerically. In Fig. [21 we show the time evo- 
lution of a displaced vortex line that was initially almost 
straight. With the parameters used in this figure, which 
was obtained by a fourth-order Runger-Kutta method, 
the vortex line is seen to collapse into three bright soli- 
tons. 

Discussion. — The in situ imagining of the three- 
dimensional vortex line is difficult, but not impossi- 
ble [3 03. However, in many experiments the mag- 
netic trap potential is turned off and the condensate is 
allowed to expand in order to facilitate the imaging of the 
Bose-Einstein condensate. In an earlier paper we 
showed that during the expansion the structure of the 
three-dimensional vortex line in a good approximation 
only changes by a scale factor and therefore vortex-line 
solitons are observable even after the expansion of the 
Bose-Einstein condensate. 

This offers the opportunity to study experimentally a 
variety of other interesting phenomena. One intriguing 
possibility is to study a Bose-Einstein condensate with 
more than one vortex line. In this manner it is possible 
to study how solitons in different vortex lines interact. 
The possibility of breather solutions is also an exiting 
prospect |20| . 
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